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Abstract 

A general multiband formulation of linear and non-linear optical response functions for realistic 
models of correlated crystals is presented. Dipole forbidden d-d optical transitions originate from 
the vertex functions, which we consider assuming locality of irreducible four-leg vertex. The uni- 
fied formulation for second- and third-order response functions in terms of the three-leg vertex is 
suitable for practical calculations in solids. We illustrate the general approach by consideration 
of intraatomic spin-flip contributions, with the energy of 2J, where J is a Hund exchange, in the 
simplest two-orbital model. 
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From the physical point of view the most of natural mineral dyes are the Mott or charge- 
transfer insulators and their colors are determined essentially by correlation effects^. If a 
material has no energy gap or its value is smaller than the energy of visual light, it will 
be non-transparent, a black one or with metallic shine. In the case of broad-gap materials 
the absorbtion of visual light is determined by impurities (ruby, that is, AI2O3 doped with 
Cr 3+ ions is a prototype example-) or by optically (dipole) forbidden d-d transitions between 
different terms and multiplets belonging to the same d n - configurations of transition metal 
ions for a pure system. The latter processes are responsible for a green color of NiO^ and 
blue color of most of the divalent copper insulating compounds^. These systems are usually 
colored and transparent and the transparency itself is a manifestation of dipole-forbidden 
character of relevant optical transitions. 

Up to now the optical properties of Mott or charge transfer insulators are considered 
within the framework of cluster approaches 4 ' 5 . The present paper develops a general trans- 
lationally invariant formalism to treat the d-d transitions in strongly correlated crystals. It 
is commonly accepted now that the standard LDA(GGA) approach is insufficient to describe 
the electronic structure of the Mott insulators 6 . To have more adequate picture of single- 
electron spectra various approaches have been applied to the problem such as LDA+U-, self- 
interaction corrections^ and the GW-scheme^. However, all these approaches do not provide 
the correct atomic limit and in particular do not take into account the term and multiplet 
structure, which is crucial for optics. This problem can be solved within the LDA+DMFT 
(dynamical mean-field theory)^ 1 ^ and the Hubbard-I approximation^ 1 ^. There were sev- 
eral attempts to calculate the optical properties within the LDA+DMFT using the Kubo 
formula for optical conductivity^— . However, in all these calculations the vertex contribu- 
tions were not taken into account and the two-particle Green functions were calculated as 
a convolution of two single-particle Green functions (for review of the LDA+DMFT see^). 
The latter contains only transitions related to the promotion of d-electrons to the p-band 
with the change of transition-metal configurations from d n to d n±1 and thus this approach 
is not sufficient to explain why NiO is green. 

In the DMFT approach the self-energy is local which leads to a cancellation of vertex 
corrections in the single-band Hubbard model 17 . However, this is not the case for a generic 
multiband situation, similar to a treatment of optical properties of disordered alloys in the 
coherent-potential approximation (CPA)^. Here we present the corresponding formulation 
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for linear and non-linear optical response functions. 

We start with the general expression for linear optical conductivity 

e 2 

cr ab (iu) = — T 2 y j } j (4 \v a \ 1) (2 |u 6 | 3) xvmii^, iu) (1) 
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where X1234 is a generalized susceptibility in a quasiatomic basis set |1) = \iLa), where i, L, a 
label sites, orbital quantum numbers and spin projections, respectively, v a is the electron 
velocity operator (a = x,y,z), and we use the Matsubara Green functions; at the end of 
calculations the analytical continuation iu — > u + iO should be performed^. 

The two-particle Green function x is expressed in terms of the single-particle Green 
function G and the irreducible vertex function 

Xi234(^, iv', iu) = -G 12 (iv')G u (iv)5 v y +U} - (2) 
-T /;/; G 15 (iu')G 84 (iu)T 5678 (iu, iu", iu)x%2m{iv" ', iv', iu) (3) 
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which can be written in the matrix form in the fermionic Matsubara frequencies (iu, iu') and 
the pairs of electron quantum numbers (14, 23) as 

X = Xo + Xof X (4) 

where xo — —G*G. Within the DMFT approximation the self-energy H(iu) is local, that is, 
diagonal in site indices and k-independent in the momentum representation^. In addition 
we will assume a locality of the irreducible vertex function T; this is the only approximation 
we add. Than, it can be obtained from the local version of Eq. 

r = xTolo - Xioc ( 5 ) 

where Xioc,o, Xioc are matrices in the Matsubara frequencies and pairs of orbital and spin 
indices, all site indices are supposed to be the same^. Both single-particle and two-particle 
on-site Green functions can be found numerically exactly using full diagonalization scheme 2 ^ 
or continuous-time Quantum Monte Carlo method 2 ^. 

Thus, the optical conductivity flTJ can be expressed in the following form: 

e 1 

a ab (iu) = T VVV(4k \v e a ff (iu, iu) | lk> G 12 (k,iu) (2k \v b \ 3k) G u (k,iu + iu) (6) 
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Figure 1: (Color online) Diagrammatic representation of Eq.Q; thick lines are the exact Green 
functions, dot is the bare matrix element (lk|v a |2k) and dot with triangle corresponds to 
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where bl234) are orbital and spin indices only and the effective matrix element (three-leg 
vertex) satisfies the equation 

(4k\v e J f (iv,iu)\lk) = (4k|v a |lk) -Tj2J2J2( 3k '\ v * f ( iu '^^\ 2k ') ^ 

v> k' 2356 

G 2 5(k', iv')G 63 (k', iv' + iu)T 5146 {iv' , iv, iu) 

Diagrammatically Eq.([6]) is shown in FigfTJ 

The effective matrix element (lk \v e J^{iv, iu) | 2k) is convenient since its use allows us 
to present the linear and non-linear response functions in the unified form. For example, 
the non-linear optical susceptibility describing a second harmonic generation can be exactly 
represented as a sum of diagrams shown in Figsj2^i and [2)3. 

The corresponding analytical expression for Figj2^i reads 

XabcilU, iu, 2iu) = Y, J2 ( 6k iuJ ) I lk ) G M k >™) ( 8 ) 

U v k 123456 

2k {iv + iu, iu) | 3k^ G 34 (k,zz/ + iu) (4k \v e J^ {iv — iu, 2iu) | 5k) G 56 {k,iv — iu) 

The calculations of magnetic susceptibility in the one-band Hubbard model^ shows that 
the contributions of the six-leg vertex are small; one can hope that this is also the case for 
the second harmonic generations and thus Eq.flS]) will be sufficient for real calculations. 

We proceed with the multiband Hubbard model with the Hamiltonian 

H = e i c ima C imo~ ^ ] ^mm' C ima C jm'a- + ^ ^ ] ^ rn\m^mz'm4^tm\a C irn^a' C ^^a' c imza (9) 

ima iy^j iacr' 

mm'cr m\rnim,zrn± 





(a) (b) 

Figure 2: (Color online) Diagrammatic representation of the non-linear susceptibility for second 
harmonic generation; thick lines are the exact Green functions, dot with triangle corresponds to 
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2k ) and the shadowed hexagon is the irreducible six-leg vertex. 



To clarify a physical meaning of vertex corrections to the response functions we discuss 
first exactly solvable model of two sites with two orbitals (i = (1,2), m = (1,2)). The 
corresponding rotationally invariant interaction matrix is parametrized by the Hubbard 
energy U and the Hund exchange parameter J as^ 

U rni m2mim,2 U 2</ 

Knim2tn2mi ^ (1^) 

(mi m-z). The dimension of the Hilbert space is equal to 2 8 so x an d G can be easily 
found by exact diagonalization. The results for lmx/oj oc Kea(u) for real (not Matsubara) 
frequencies are shown in Figj3]and FigJH 

As one can see from the density of states the single-particle transitions which manifest 
themselves in xo starts at the frequency ou > 1 corresponding to the distance between the 
highest occupied and the lowest unoccupied orbitals. For small enough hopping (Figj3^) 
there is a lot of peaks at smaller energies which originate from the poles of the vertex 
functions and represent the optically forbidden transitions in our toy model. In particular, a 
transition with u = 2 J = 0.4 is clearly visible for intraband intersite transitions (11) — > (21). 
The transition corresponds to local intraatomic spin flip processes. It is visible in optical 
(summed up over spins) susceptibility, since the Coulomb interaction matrix couples spin-up 
and spin-down states. An interesting and unexpected result of the toy model is that for a 
moderate hopping (Figj3b) all these local term effects disappear. 
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Figure 3: (Color online) Components of generalized susceptibility Imximjm'^m.jm'/^ summed up 
over (a, a') for two-site two-band model with E\ = e% = 0, U = 1, J = 0.2, ti2 = £21 = 0.05; (a): 
tn = ^22 = 0.2; (b) in = ^22 = 0.5. Solid red curve corresponds to intrasite interband transitions 
im = (11), jm' = (12), dashed blue curve - to intersite intraband transitions im = (11), jm' = (21), 
and the dotted green curve - to intersite interband im = (11), jm' = (22). We use the spectral 
representation with the Imu; = 0.05 and the temperature is T = 0.05. Energies are in the units of 
U. 

Now we present the main part of our work related with the two-band lattice model with 
the use of local approximation for vertex f as describes above. The calculations have been 
done for the square lattice in the nearest neighbour approximation. The local vertex function 
has been obtained from Eq.([5]) by exact diagonalization calculations of x and G. If the exact 
diagonalization would be performed for intraatomic Coulomb interaction Hamiltonian only, 
this would correspond to the Hubbard I approximation 11 . To go beyond this, in spirit of the 
DMFT 17 , we have added one more orbital to the bath. 

Despite simplicity of the model our calculations turned out to be rather time and mem- 
ory consuming due to the inversion of x _m atrices depending of two Matsubara frequencies 
[iv,iv'\ To reach convergence, we had to use about one hundred frequencies and 20 x 20 
k-points. The computational results are shown in Figj5] One can see that £0 has only one 
pronounced peak at u = U corresponding to the transition from the lower to upper Hubbard 
bands. At the same time, x has an additional peak within the gap of the single-particle 
excitation spectrum at the maximum around u = 2 J = 0.4. This maximum originates from 
the spin-multiplet structure of the d 2 configuration. 
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Figure 4: (Color online) Density of single-particle states for the same parameters as in Fig. 3 
Energies are in the units of U. 
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Figure 5: (Color online) Intrasite, interband components of Imx/w oc Re<r(u;) (solid red curve) and 
Imxo/w oc Recro(w) (dashed blue curve) for the square lattice in the nearest-neighbour approxima- 
tion, with the same parameters as in Fig. 3a. Energies are in the units of U. 

To conclude, we present a general formalism which allows to consider term and multiplet 
effects on linear and non-linear optical properties of multiband strongly correlated systems. 
For a two-band model the computational results are quite reasonable and one can hope that 
the scheme can be applied to the first-principle investigations of realistic systems in the 
spirit of the LDA+DMFT. The approach can be apply for other response functions such as 
magnetic susceptibilities and STM-spectroscopy where the transition with energy u = 2J 
have been recently observed for Mn-chains on Pt surface^. 
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